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Abstract. The minimal nearest-neighbor tight-binding model with the Peierls substitution is em- 
ployed to describe the electronic structure of Bemal-stacked graphite subject to tilted magnetic 
fields. We show that while the presence of the in-plane component of the magnetic field has a neg- 
ligible effect on the Landau level structure at the K point of the graphite Brillouin zone, at the H 
point it leads to the experimentally observable splitting of Landau levels which grows approximately 
linearly with the in-plane field intensity. 
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INTRODUCTION 

Recently, graphene monolayers have attracted much attention motivated by their un- 
usual two-dimensional (2D) Dirac energy spectrum of electrons. In graphite, which is a 
graphene multilayer composed of weakly coupled graphene sheets, the interlayer inter- 
action converts the 2D electron energy spectrum of graphen into the three-dimensional 
(3D) spectrum of graphite. The application of the tilted magnetic field is a typical method 
to distinguish between 2D and 3D electron systems, as in 3D systems the orbital effect 
of the in-plane magnetic field should be observable. This problem has been touched in 
two recent theoretical studies of the graphene multilayer energy spectrum in magnetic 
fields parallel to the layers [1], and of the Landau levels (LLs) of bilayer graphene in 
magnetic fields of arbitrary orientation [2]. Both papers conclude that the very strong 
in-plane field components are necessary to induce observable effects on the electronic 
structure. Here we show that at the H point of the hexagonal Brillouin zone of graphite 
the application of the tilted magnetic field leads to experimentally observable splitting 
of LLs. 



MODEL AND RESULTS 



Bulk graphite is composed of periodically repeated graphene bilayers formed by two 
nonequivalent Bemal-stacked graphene sheets. There are two sublattices, A and B, on 
each sheet and, therefore, four atoms in a unit cell. The distance between the nearest 
atoms A and B in a single layer is 1.42 A, the interlayer distance between nearest atoms 
Ais J = 3.35A. 



To describe the graphite band structure, we employ the minimal nearest-neighbor 
tight-binding model introduced in Ref. [3] and successfully applied in Refs. [1, 2, 4]. 
The wave functions are expressed via four orthogonal components i//^, i^J, vf+i' '^f+i' 
which are, in zero magnetic field, Bloch sums of atomic wave functions over the lattice 
sites of sublattices A and B in individual layers. The tight-binding Hamiltonian ^ 
includes only the intralayer interaction — 3.16eV between nearest atoms A and B 
in the plane, and the interlay er interaction ? = 0.38eV between nearest atoms A out 
of plane. The continuum approximation is used in the vicinity of the H — K — H axis 
of the graphite Brillouin zone, for small k measured from the axis. Then the electron 
wavelength is larger than the distance between atoms, and the non-zero matrix elements 
of J^f can be written as 

jif^B = hvF{kx + iky), M'^^ = hvF{kx - iky), = t. (1) 

The Fermi velocity, vp, is defined by hvp = ^/^ajQ/l and will be used as an intralayer 
parameter instead of 70 in the subsequent consideration. 

The effect of the magnetic field B — {Q,By,B^ is conveniently introduced by Peierls 
substitution into the zero-field Hamiltonian. If we choose the vector potential in the form 
A = {Byz - B^y, 0, 0) we get 

hkx — >■ fikx — I e\Bzy -\- \ e\Byjd. (2) 

When solving the Schrodinger equation for the above four wave functions, we can 
make use of the structure of J^f, which allows to express the function xj/? via the function 

X}/j from the same layer, and leaves us with two „interlayer" equations for and V^^+i- 
Note, that graphite is no longer periodic in the z-direction (parallel to the c-axis) but 
becomes periodic in the direction of the tilted magnetic field. Here we apply an approach 
developed in [5]. The new periodicity implies that i/zj^ can be written as 

Xjff = e''^'j(^t (7] +;77d) , V^+i = e^'^'^^J+'^i (t] + {j+l)r]d) , (3) 

—n/2 < k^d < n/2. Here {rj + jrjd) and 02 + Jld) can be associated with cy- 
clotron orbits in two layers. The dimensionless variable rj is defined by r] = {y — yo)/£, 
where _yo is the cyclotron orbit center, jo = -^jkx, and 4 is the magnetic length, = 
h/{\e\Bz). The small parameter 

By d 

means the shift of the cyclotron orbit center in the j-layer due to the in-plane component 
of the magnetic field. By. Introducing the shift operator by 

0(r] + r]rf) = e''^''''0(7]), K=-id/dn, (5) 
and employing JC-representation, two „interlayer" equations can be given the form 

(l>2 +2tE cos {Krid + kzd)(^f ^0, (6) 



+ K^-l] -E 



+2tE COS (Kria + k^d) (1)2 =0, (7) 



where ^ stands for 7}i\e\v\B^. Eqs. (6,7) represent the main resuh of this work. 

It is evident that 0^ and 02 cannot be too far from the eigenfunctions of the harmonic 
oscillator, ^„(k"). If we write 
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(8) 



n'=0 



n'=0 



we amve to 



[^(n+l)-£2]^2,n + 2£ £ VAiy =0, 

OO 

[^n - £2] Ai,„ + 2£ £ 7;,„,A2,„/ = 0, 



n'=0 



where 
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(10) 

(11) 



In the perpendicular magnetic field = ^i^d r„ „/ = 2 cos {kzd)5„^„i. This corre- 
sponds to the previously discussed case [3]. At the K point, k^d = 0, we get the energy 
spectrum of an effective bilayer. At the H point, k^d = ±;r/2, the coupling between 
layers disappears, and we obtain the LLs corresponding to graphene Dirac fermions, 
namely E^^ = ±.\fWn for the first layer, and E^^^ = ±a/^(« + 1) for the second layer, 
n = 0,l,2,--. 

Let us discuss the magnetic field of an arbitrary direction. It is obvious that the 
influence of the small parameter 7]^ is negligible at the point K when k^d = and 
cos [kt] J — k^d) fa 1. This case corresponds to the bilayer subject to the tilted field 
discussed in Ref. [2], with the result that corrections induced by By are negligible. 

Let us consider the influence of By at the H point, when cos{kzd) = and sm{kzd) — 
±1. Then, unlike the case of the perpendicular field, the interlayer interaction is not 
reduced to zero, but remains finite. Note, that the small perturbation r/^ couples the un- 
perturbed states j/ij) and |(n+ l)f ), which belong to the same unperturbed eigenvalues 
±-y/^(n+ 1). Consequently, the perturbation approach suitable to describe the coupling 
of these degenerated states must be applied, which yields equations 



[^{n +l)-E^] A2,n + 2ETn,n+lAi,n+l = 0, 

[^{n+l) -^2] Ai,„+i +2Er„+i,„A2,„ = 0, 



where 



n + l 




+ . 



(12) 
(13) 

(14) 



The secular equation derived from Eqs. (12, 13) reads 



[j3§{n+l)-Ef-4E^T^^„^,=0, 



(15) 




FIGURE 1. LLs of graphite at the H point as a function of the perpendicular component of the magnetic 
field, B^. a) LLs, Zij^^^ and E^^, in the perpendicular magnetic field {(p = 0")- b) LLs, E^'"^, splitted by 
tilting the magnetic field ((p — 20"). 



and from here we get the four eigenenergies 

£„^;t = ±r„,„+i ± y/^H^iy+T^, n = 0,l,2,---. (16) 

The eigenenergies, and , originated from E^^ remain the same as in the per- 
pendicular magnetic field. In that case the degeneracy is not removed. Lifting of LL 
degeneracy by the tilted magnetic field is shown in Fig. 1. The LL splitting is of the 
order of several meV, and it grows with the tilt angle, i.e., with By. 

To conclude, we have obtained the LL structure of graphite at the H point of the 
Brillouine zone in the tilted magnetic field configuration. This effect is experimentally 
observable. 
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